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ABSTRACT 

We analyze the question of screening versus confinement in bosonized massless 
QCD in two dimensions. We deduce the screening behavior of massless SU{Nc) 
QCD with flavored fundamental fermions and fermions in the adjoint representa¬ 
tion. This is done by computing the potential between external quarks as well 
as by bosonizing also the external sources and analyzing the states of the com¬ 
bined system. We write down novel “non-abelian Schwinger like” solutions of the 
equations of motion, compute their masses and argue that an exchange of massive 
modes of this type is associated with the screening mechanism. 
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1. Introduction 


The question of confinement versus screening in four dimensional (4D) non- 
abelian gauge thoeries is one of the major problems of high-energy theory. Two 
dimensional (2D) gauge theories may serve as a laboratory in exploring that prob¬ 
lem. Just as in 4D, also in 2D one can use the potential between two heavy external 
charges, the expectation value of a Wilson loop and the structure of the bound state 
spectrum as probes of conhnement. It is believed that a potenial growing linearly 
at large separation distances, an area law behaviour of the Wilson loop and a spec¬ 
trum which is independent of the number of colors W indicate that the system is 
in a conhning state. Note that in some cases, like with Higgs in the fundamental, 
screening and confinement are in one phase by “complementarity” 

In a paper by D.Gross et al. it was argued that there is a screening effect 
between heavy external charges induced by massless dynamical fermions even if 
the latter are in a representation which has zero “iVc-ality”, namely, vanishing 
center ~ 0). It was further shown that conhnement is restored as soon as 

the dynamical fermions get some non-trivial mass. In that paper both the nature 
of the potential and the Wilson loop were determined in the abelian theory and 
in several non-abelian cases. In one case the group was SU{Nc = 2) with the 
dynamical fermions in 3, and the other case was SU{3) with 8 fermions. In the 
latter it was shown in fact that the spinor 8 of SO(8) are screening. 

The potential between the external quarks can be extracted using several dif¬ 
ferent methods: (i) Deriving the effective Lagrangian (integrating out the fermions) 
and then extracting the potential using the static gauge conhgurations that solve 
the corresponding equations of motion, (ii) Using the gauge configuration that 
solves the equations of motion of the bosonized action . (iii) Eliminating the 
gauge helds from the bosonized gauged action, solving for static solutions of the 
currents and deducing the potential as the difference between the Hamiltionian of 
the systems with and without the external sources. 

In the present paper we use method (ii) to prove the statements of screening 
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for massless dynamical fermions. 

Another approach to determine if the system conhnes is based on bosoniz- 
ing also the heavy external charges!^^ Conhnement manifests itself in this double 
bosonized model by the absence of soliton solutions that correspond to unbounded 
quarks. In case that there are quark hnite-energy static solutions, one may con¬ 
clude that the system is non-conhning. We use also this type of analysis for both 
abelian and non-abelian gauge theories. 

The screening mechanism in the massless Schwinger model could be attributed 
to the exchange of the emerged massive photon, which is the only state in the 
exact spectrum. The non-abelian counterpart is clearly much more complicated 
and seems to be a non-integrable model However, by introducing flavor degrees 
of freedom one can pass, in the limit of large number of flavours Nf (with hnite 
Nc), to a domain where the non abelian theory resembes a collection of — 1 
abelian theories. In that limit the spectrum includes — 1 massive modes of the 
type that exist in the Schwinger model . One can then draw an intuitive picture 
of screening due to those modes in a similar manner to the one in the abelian 
theory. As a matter of fact, it is only in the large Nf limit, that one can justify 
relating solutions of the equations of motion and physical states and deducing 
conclusions about the sctructure of the spectrum. One might hnd in the “massive 
gauge states” an indication of the “non-conhning” structure of the spectrum. The 
reason for that is, that even though they are gauge invariant states, they are in the 
adjoint representation of a “global color symmetry” and not singlets of that group. 
These states had already been pointed out in an earlier workj*^'based on a BRST 
analysis and a special parametrization of the gauge conhgurations. However, in 
that paper we were not able to rigorusly show that indeed they were part of the 
BRST cohomology. Note, that even if they are not in space of physical states, the 
massive states could nevetheless be responsible for the screening potential. 

When passing from a screening picture at large Nf and hnite Aj, to the domain 
of a small number of havors, one can anticipate two types of scenarios: (i) A smooth 
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transition where the screening behavior persists all the way down to Nf = 1; (ii) 
A phase transition at a certain valne of Nf and a confining natnre below it. One 
may argue that the massive modes of large Nf are an artifact of the abelianization 
of the theory. To check that possibility we have searched for non-abelian solutions 
of the equations of motions. Indeed, we found new non-abelian solutions that are 
also massive and are associated with the gauge helds, namely, are in the adjoint 
of the global color group. We thus conclude that this nature does not stem from 
the abelianization of the large Nf limit, and hence may present certain evidence 
in favour of option (i). A different “patch” of the space of physical states for hnite 
Nc and hnite Nf, that of the low lying baryonic states, was determined in the 
semiclassical domain in [7]. Since those baryons are all color singlets this might 
seem as a contradition to the screening nature of the spectrum. In fact there is 
no contradiction since the baryons were discovered only for massive quarks and 
not for massless ones. As was shown in [2] turning on a mass term for the quarks 
changes the picture dramatically into a conhnig one. 

’t Hooft solved the spectrum of QCD 2 in the large N(, limit. In that analysis the 
quarks were havorless and in the color fundamental representation. This procedure 
was recently also applied for adjoint fermions. Both in the original work as well as 
in those of ref. [9] there is no trace of the massive modes that our work analyzes. 
Differently stated the large W approach reveals a conhnig spectrum both for the 
case of massless and massive quarks. We believe that this is an artifact of the 
large W limit and at hnite W the spectrum of a theory with/without quark mass 
behaves like a conhning/screening spectrum respectively. 

The paper is organized in the following way. In section 2 we review the rules 
of bosonization of two dimensional QCD with both massive and massless fermions 
which transform in the fundamental or adjoint representations. The equations of 
motions of bosonized QCD 2 in the presence of external currents are derived and 
discussed in section 3. Non-abelian solutions of the equations as well as some in¬ 
teresting abelian ones are presented in section 4 for the model without external 
sources, and in section 5 when the latter are turned on. The energy-momentum 
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tensor and the spectrum that corresponds to the solutions of sections 4 and 5 are 
derived in section 6. Section 7 is devoted to the analsis of the system by bosonizing 
also the external currents. In section 8 we summarize the results of the present 
work, state our understanding of the nature of QCD2 in the different regimes 
and raise some further intersting open questions. In appendix A we hnd some 
“truly non-abelian” SU{Nc = 2) solutions of the equations of motion of bosonized 
QCD2. Appendix B is devoted to the derivation of the energy momentum tensor 
of a bosonizied fermion coupled to an abelian gauge held and its non-abelian gen¬ 
eralization. In appendix C we show that minimal energy solutions of the massive 
WZW model are necessarly in a diagonal form. 

2 . Review of Bosonization in QCD2 

We start with reviewing the bosoization formulations of QCD2 with fermions 
in the fundamental and adjoint representations. 

Dirac fermions in the fundamental representation. 

Multihavor massive QCD2 with fermions in the fundamental representation 

[7] 

was shown to be described by the following action 

Sqcd 2 =Si{u) - [ cPzTr{iu~^duA + iudu~^A -\- Au~^Au — AA) 

( 1 ) 

j if 2 '■ Trclu.+ «-'] : j d^zTrjjlF^] 

where u G U{Nf x Nc)] S^u) is a level k WZW model; 

5'fc(w) = [ (fxTr{duud'^u~^) + [ d^ye'^^^Tr{u~^diu){u~^dju){u~^dku) 

OTT J 127r J 

B 

( 2 ) 

A and A take their values in the algebra of = SU {Nq)', F = dA — dA + i[A, A]; 
w? equals niquC, where /i is the normal ordering mass and C = with 7 
Euler’s constant. Notice that the space-time has a Minkowski signature, and we 
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use a notation in which the light-cone components of a vector are denoted by 
B = i?_|_ and B = i?_. 


The action for massless fermions can be simplihed using the following parametriza- 

T h c QTTf AT \ n c QTTf AT A / c _ UjN^Nf) _ 


tion u = ghle V where h e SU{N,), g e SU{Nf) and I G xU(4)]iv.xt/^(i) • 

The action in terms of these variable takes the form^^' 


N f 

S = SNf{h) - - / (fzTr{ih~^dhA + ihdh~^A + Ah~^Ah — AA) 

2tt J 

+ Snc{ 9 ) + ^ j (fz[d(l)d(l)] ( 3 ) 

+ 4 / dhTrH[F^] 

[ 3 ] 

Notice that the action is independent of /. Recall that to discuss the quark soliton 
structure we need the massive u G U{Nf x Nc) description!"^' 

Majorana fermions in the adjoint representation. 

A non-abelian bosonization of Majorana fermions in the adjoint representa¬ 
tion'^”'can be expressed in terms of S{had) where had are {N^ — 1 ) x (N^ — 1 ) 
matrices, so that the action for the corresponding QCD 2 now reads, in the mas¬ 
sive case, 

S = ^S{had) - ^ j d^zTr{ih~^dhadA + iudh~^A + Ah~^Ahad - AA) 

+ 4 / d^zTrniF^] 

J 

The factor ^ in front of the S{had) term comes from the reality nature of the 
Majorana fermions. It is straightforward to realize that the conformal anomaly of 
this model is indeed c = — 1), the associated currents have affine Lie algebra 

with an anomaly of Nc and that the conformal dimension of had is ‘Ahad = 5 (left 
and right dimensions and a total conformal dimension of 1). 
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3. Equations of motion of QCD 2 in the presence of external 

CURRENTS 

The equation of motions which follow from the variation of the action (3) with 
respect to h are given for the massless case by 

d{h~^dh) + id A + i[h~^dh, A] = 0 

( 5 ) 

d{hdh — id{hAh ^) = 0 

where a gauge A = 0 has been chosen. Notice that the second equation can be 
derived from the hrst one by multiplying it with h from the left and h~^ from the 
right. A similar result but with had replacing h follows the variation of eqn. (4) 
with respect to had- 

As was discussed in the introduction our aim is to analyze the system in the 
presence of external sources. External currents are coupled to the system by adding 
to the action (3) or (4) the following term 

h^ext — ^ J d zTvi^JextA-\~ J^xtA'). 

The variation of the combined action with respect to A and A (and then setting 
A = 0) yields for the case of Nf fundamentals the following equations of motion 

d‘^A + ac(iNfh~^dh + Jext) = 0 

( 6 ) 

ddA + [idA, A] — q;c[A^/(^^^^ ^ + hAh ^ — A) + Jext] = 0 
2 

where etc = ^. If follows from the equations of motion (5) and ( 6 ) that both the 
dynamical currents jdy = ^^h~^dh, jdy = ^^[hdh~^ — ihAh~^ + iAC\ as well as 
the external currents are covariantly conserved, which for A = 0 reads 

Djdy T djdy = 0 DJext T dJext — 0. (7) 

with D = B — i[A,]. One can eliminate the dynamical current and derive the 
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following equation for the gauge fields in terms of the external currents 


ddA + [idA, A] + a^NfA - Jext) = 0 


( 8 ) 


In fact the equation one gets in this way is the d derivative acting on the l.h.s of 
(8) equals zero. However, one can £x the residual gauge invariance A iu~^du + 
u~^Au^ with u an anti-holomorphic function du = Q (thus preserving A = 0), to 
eliminate the antiholomorphic function that should have been put in the r.h.s of 
(8). Thus A = ih~^dh and jdy = Note that this is not the 

current of the free case which is j/ree = ^^hdh~^ = —^^{dh)h~^ 


The last equation, which holds for the massless cases, (3) and m = 0 of (4), is 
universal in the sense that it is independent of the representation of the dynamical 
fermions Once mass terms are added there is an explicit dependence on the 
dynamical fermion and instead of (8) one finds that the d derivative of its left hand 
side is equal to im^Nfa^h — h~^), where m? is given in eqn. (2) and where for the 
fundamental representation we have chosen non-flavored configurations namely we 
set of eqn.(3) to 1,1,0 respectively (otherwise instead of a factor Nf, h is 


multiplied by the flavored contribution Trgle 



Studying the quantum system by analyzing the corresponding equations of 
motion is a justified approximation only provided that the classical configurations 
dominate the functional integral. Such a scenario can be achieved for the case 
of massless quarks in the fundamental color representation in the limit of a large 
number of flavours. In fact one can show that in that case ^ plays the role of /h. 
The colored sectors of the action takes the form 


S = Nf\ Si{h) - — I d^zTr{ih-^dhA) + ^ / d^zTr[dA] 


2ti 


(9) 


where Cg = Cc^jNf. 



4. Solutions of the equations without external quarks 

Let us consider first the case where the exteral sources are switched off. It is 
obvious that an “abelian” massive mode is a solution of the equations of motion 
( 8 ). Consider a conhguration of the form A = T‘^A°'{z, z) = T°'d°''°‘°A{z, f),where 
ao is a given index that takes one of the values 1 , — 1 , then the commntator 
term vanishes and A has to solve ddA + olcA = 0 with oic = Nfac- It is clear that 
there are — 1 such solutions and in fact it is easy to see that this property will 
be shared by every possible solution. This follows from the fact that the eqnation 
of motion is not invariant bnt rather covariant with respect to the “global color” 
transformation A u~^Au with a constant u. 

Let us now check whether the eqnations admit soliton solntions. For static 
conhgurations the eqnation reads 

dlA - V2[idiA, A] - 2acA = 0 (10) 

Multiplying the equation by A, taking the trace of the result and integrating over 
dx one hnds after a partial integration that f dx[rr[{dAY + 2 q:cM^] = 0 which can 
be satished only for a vanishing A. 

In the search for other possible solutions one may be instructed by the fact 

that Fzz = dA can be written in 2 d as Fzz = ^zzF and impose an ansatz for the 

solution of the form A = zC{p) where p = zf. Expanding C as a power series 

in p one hnds that the commntator term has to vanish and C is determined by 

the equation pC” + 2C' + acC = 0, where C = dpC. By a change of variable 
1 

C = p 2 lF(x) one can rewrite this eqnation as IF” + Fy/' + (1 — ^) 1 F = 0 with 
x^ = 4cicd, which is x~‘^ times a Bessel eqnation of order one, so that the solution 
for the gange held takes the form 

A = ^2^jA2V^) ( 11 ) 

yaczz 

where Aq is an arbitrary constant matrix. 
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The next task is to examine whether there are any possible solutions which 
are “non-abelian” in their nature. Consider in the special case of SU ( 2 ) the con- 
hguarion A = with a constant matrix Aq = cqTq + er + er . Plugging 

this ansatz into eqn. ( 8 ) with no external source one finds that there is a solution 
provided that 86 and 86 are constants, namely 6 = 6q + kz + kz where k, k and 
6*0 are constants. Indeed the following gauge field 


.j an — kk 


{kk - dc)dcr -ie 


[e-*V + A^f] 


Ak^ 


( 12 ) 


is a “ non-abelian solution”. The notation and the derivation are presented in 
appendix A. Setting = 0 requires that {kk — > 0. Looking into the case 

where Aq is not a constant, but with 8Aq = 0 to preserve F{Aq) = 0, one gets 
that the only solution is with a constant Aq 

The solution we have found eqn. (12) are truly non-abelian solution. The 
corresponding F is F = ike~'^^'^° {er — er)e*^'^°. Performing a gauge transforamtion 
with U = one hnds that Fjj = ik{eT — er), Ajj = Aq + Zero and Ajj = kro 

(recall that A = 0). We thus found that Ajj^Ajj and Fjj are fixed in space-time 
and no two commuting. Furthermore an abelian gauge conhguration of the form 
A = —ikz{eT — er) and A = 0 that leads to the same F is not connect to Ajj, Ajj 
by a gauge transformation. 

Using the expression of A one can easily extract jdy and jdy This will be done 
in section 6. Moreover, one can determine the non-abelian group factor h. From 
eqn.(6) it follows that h~^8h = Using the ansatz for A = it is 

easy to hnd that 

= g-Adro+e(r+r)le*ero .^3^ 


where e 



1 . 
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5. Solutions of the equations with external current 

Next we want to turn on a covariantly conserved (eqn. (7)) external current 
Jext and study the corresponding equations of motion. Abelian solutions are easily 
constructed. For instance for a pair of quark anti-quark as an external classical 
source J°‘ext = — R) — 5{xi + i?)] the abelian solution is 

A = (14) 

Inserting this expression into ^ f dxiTr[AJgxt] one hnds the usual screening po¬ 
tential 

V{r) = - e-2F^I^I)Tr[(T“°)2] (15) 

271 

Again the challange is to hnd “non-abelian” solutions where the commutator 
terms do not vanish. The SU (2) “non-abelian solution” of above is a solution also 
in case of a constant external current J = with the trivial modihcation 

/ ®c[^ + (fcfc —Oc)] 

that g ( see appendix A) is replaced by u--. Consider now an external 

current of the form J = Jo(7)ro. A solution in that case is 

A = (Jo + Jo(z))ro + + c.c] 

where dl{z) = with /o and g related to k and k as given in eqn. (12). In 

the case of light-front quatization with ^ playing the role of the space coordinates, 
JextA) stands for a general “static” current. In particular a current density that 
corresponds to a quark anti-quark pair takes in this framework the form = 
— R) — d{z + i7)] and the corresponding solution has e{z — R) and 
e{z + R) factors in 6. The corresponding potential is a constant thus non-cohning. 
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6. The energy-momentum tensor and the spectrum 


Next we want to analyze the spectrnm of physical states that correspond (at 
least in the large Nf limit) to solutions of the equations of motion. Recall that those 
states transform in the adjoint representation of the global color transformations. 

First we have to compute the energy momentum tensor T = Tgz,T = T^ziT^z 
that corresponds to the action (3). Only the colored part of the energy momentum 
tensor is relevant to our discussion. From appendix B 


T = -- "" -- : TrluyJdy] : T.z = ^Tr[{dAf 


(16) 


Nf + Nc SttcTc 

where the currents of the dynamical quarks which were dehned below eqn.(6) are 

(17) 


iNf _i 

jdy =-r^h ^dh = 


Jdy 


271 

1 


27Ta, 


2'Ka, 


■{ddA -I- i[dA, A]) = 


iN 


Nf^ 


277 


^h-^dh = -'-^A 

277 


To proceed and compute the masses of the physical states one has to choose 
a quantization scheme. It is natural in the light cone gauge to use a light front 
quantization. In that scheme we take 2 ; to denote the space coordinate. In appendix 
B we express the momentum coponenets P and P as integrals over T and Tzz- The 
masses of the states are given by the eigenvalues of = PP. 

To set the proper normalization of the helds let us consider hrst the abelian 
solutions for A. In that case the operator A can be written as 


Ts = r^y +C.C] 

where is a matrix in the Cartan sub-algebra, kk = ac, the creation and annihi¬ 
lation operators obey the commutation relation [a{k), A (k)] = 6{k — k) and J\f{k) is 
a normalization factor. Inserting this form of R it is a straightforward calculation 
to get J\f{k) = 2 ^/ 7 ^ so that on the states \k,k > is equal as expected to die- 
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One can instead assume a finite system of size L in z direction. In that case 
if one uses a normalization where kL is an integer n times 2tt, and that Aah = 
-^ksinO (in hrst quantized version). One then hnds that P = k and P = k, so 
that again = cuc- 


In case of the SU{Nc = 2) non-abelian solutions the superposition principle 
does not apply and there is no room for Fourier expansion of helds. We thus 
invoke the second quantization method of above, with a finite size system. The 
expectation values of the energy momentum tensor in those states are determined 
by substituting the expression for the gauge conhguration A into eqn.(16). This 
leads to 


<T > = 


<Tzz> = 


8n Nf + Nc 

I^L^kk){i-^) 

167r^ kk’ 


etc ^ 

w 


( 18 ) 


If we use the same quantization scheme as for the abelian case we get that 


P 


1 

Stt Nj + Nc 


- -kk 
{Lk)k^{l 
etc 





where L is the size of the system. Taking again the normalization kL = 2nn one 
hnds that the non-abelian state is characterized by masses 




rP Nj 

32 Nf + Nc 


dc 


dcN2 

kk’ 


The discussion above was all for Oq = 0. Thus we require that {kk — dc) > 0 (the 
case of zero kk — oic corresponds to vanishing A). We get an M starting from 
M = 0 and growing up linearly in kk for kk » dg. Note also that our solution is 
singular for Lc = 0. 


13 



7. Bosonized external currents 


Another approach to the coupling of the dynamical fermions to external cur¬ 
rents is to bosonize the “external” currents. Let us briefly summarize hrst the 
abelian case. Consider external fermions of mass M and charge qe described by 
the real scalar hied <h together with the dynamical fermions of unit charge e and 
mass m associated with the scalar (j). The Lagrangian of the combined system after 
integrating out the gauge helds is given by 

C =^((9^0(9^^) -f mTi[cos{2y/^(j)) — 1] + -I- ME[cos(2v^*h) — 1] 

(19) 

Let us look for static solutions of the corresponding equations of motion with hnite 
energy. Take, without loss of generality, 0(—cx)) = <h(—cx)) = 0. From the M term 
we get <h(cx) = ^/^iN with N integer. For m 7^ 0 we also get (f){oo) = Now 

from the term, n + qN = 0. Thus, for instance for = 1 hnite energy solutions 
occur only for q = —n. In the massless case we have only <h(cx) = ^/^iN and then a 
hnite energy solution for = 1 is if 4>{oo) = —y/nq. So, when q = —n, the system 
is always in the screening phase, whereas when q 7 ^ —n it is in the conhnement 
phase for m 7 ^ 0 and screening for m = 0. 

Proceeding now to the QCD case one can consider several diherent possibilities 

dynamical fermions that can 
be either massless or massive. Obviously, for the external source one would assign 
a mass which should then be taken to inhnity. The system of dynamical adjoint 
fermions and external fundamental quarks can be described by an action which is 
the sum of (4) and (1). 

Integrating over the gauge degrees of freedom one is left with the terms in (1) 
and (4) that do not include coupling to gauge held together with a current-current 
non-local interaction term. For the interesting case of dynamical quarks in the 
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adjoint and external in fundamental we get for the interaction term 



( 20 ) 


where u is dehned in (1), Tp are the SU{Nc) generators expressed as {NcNf) x 
{NcNf) matrices in the fundamental representation of U{Nc x Nf) and the 
(N^ — 1) X (N^ — 1) matrices in the SU{N(.) adjoint representation. The other 
cases can be treated similarly. For simplicity we discuss from here on the case of a 
single flavor. 

Let us hrst consider the case of external adjoint quarks Uext{x) G SU{Nc) x 
Ub{^), which can be represented as 


Uext 



1 




( 21 ) 


V 


l/ 


(<h is not normalized canonically here). The reason that we take a diagonal ansatz is 
that it corresponds, as we argue in appendix C, to a minimal energy conhguration. 
Ansatz (21) corresponds to QlxtQext^ namely, to an external adjoint state. We 
expect this state to be screened by the adjoint dynamical fermions. With this 
ansatz ^{ih~^^dhext) takes the form 


$ 




\ 


( 22 ) 


It is thus clear that only Tp contributes to the trace in (20). To show the dynamical 
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configuration that screens, take 


log had 


. 0-00 
0 0 0 
0 0 0 


\ 


V 



with the matrix that contributes to the Tr in ( 20 ), 


Tadj 



1 

0 

0 


0 

0 

0 


0 


\ 


V 


0 ’ 


(23) 


(24) 


which corresponds to the generator of rotation in direction 3 for the sub 0(3) of 
hrst three indices, thus obtaining the a term proportional to (<h + 0 )^ emerging 
from (20). The mass terms for u^xt and h are now proportional to (1 — cos<h) and 
(1 — COS0) respectively. A boundary condition <h(cx)) = 27r can be cancelled in the 
interaction term by the boundary condition 0(cxd) = —27r. 


Let us examine now the case of a single external quark 


'^ext 


( 


e 


—i>I> 


1 


V 


1 


\ 

/ 


(25) 


Its contribution to the interaction term is ^{Vi^ + ^^dyn”)"^ where r]i = 

, 0 and ^’’dyn” is the part of the dynamical quarks. If again we take a con- 
y 2i(i+l) 

hguration of the dynamical quarks based on a single scalar like in (23) we get 
altogether an term of the form ( 5 $ + 0)^. Now if <h(cx)) = 27r one cannot hnd a 
hnite energy solution since from the mass term 0 (oo) = 27rn, and thus there is no 
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way to cancell the interaction term. If, however, we consider massless dynamical 
fermions there is no constraint on (p{oo) so it can be taken to be equal —tt, and thus 
again a screening situation is achieved. This argument should be supplemented by 
showing that one cannot hnd another configuration besides (23) that may cancel 
the ? 7 j<h term in (? 7 j<h + “dyn”)^, for SU{Nc) with Nc > 3. 

8. Discussion 

In the present paper, using bosonization techniques, we have presented further 
evidence for the non-abelian screening of external charges by dynamical massless 
fermions. We have shown it explicitly for dynamical fermions in the adjoint and 
Nf fundamental representations. In fact, the latter case implies that a WZW 
model coupled to non-abelian gauge helds is a screening model for any level of 
the affine Lie algebra. The fact that there is no relation between the charges 
of the screening dynamical fermions and those of the external sources may seem 
unintuitive. However, one can understand this phenomena in a simple way if one 
realizes that the interaction between the external charges involves an exchange of 
a massive mode which is an outcome of the dynamics of bare massless quarks. The 
main outcome of the present paper is the observation that indeed such massive 
modes manifest themselves in the form of soutions of the equations of motion. 
This is well known for QED 2 , and has been emphasized more recently for QCD 2 ^^^ 

A natural question to ask is whether there are consequences of the screening 
behaviour in the spectrum of the theory. A simple minded intuition of the differ¬ 
ence between a conhning and a screening spectrum can be derived from quantum 
mechanics. A potential of the form (15) leads to a spectrum of bound states with 
energy smaller than the asymptotic value of the potential. A linear potential, 
on the other hand, can accomodate an inhnite spectrum of bound states with no 
limitations on their energies. Practically, of course, higher energy states will be 
unstable. It is obvious that the massless Schwinger model which has a single state 
in its spectrum falls into the former class. It looks plausible that the spectrum in 
the non-abelian case is also limited. 
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Another way to distinguish between conhning and non-conhning spectrum is 
the dependence of multiplicity of the physical states on the number of colors. In 
a conhning spectrum one hnds only color singlets and their multiplicities do not 
depend directly on N^. The massive states discussed in the present work, both the 
abelian and the non-abelian solutions, admit a degeneracy of — 1, or stated 
differently those states are in the adjoint of the “global color symmetry”. 

We have used an argument that in the large Nf limit the classical solutions 
of the equations of motion dominate the functional integral. However, it is not 
obvious that the corresponding states are physical. In a previous paper we have 
used a special formulation of QCD 2 in terms oi A = if~^df, A = ifdf~^ with 
f{z,z), f{z,z) E [SU{Nc)Y the complexihcation of SU{Nc). Implementing BRST 
techniques we have solved for the physical states of the abelian analog. In the non- 
abelian case one hnds massive color singlet states with a mass of Cc-y/^ which 
are the analogs of the states discussed in the present paper. Unlike the massless 
partners of these states, for the massive ones we were not able to show that they are 
not physical states. We still do not have a dehnite answer to this question; however, 
the fact that in the abelian theory of the large Nf limit they are physical states 
supports the conjecture that they are physical also for hnte Nf. Furthermore, it 
seems that no matter whether they are in the sub-space of physical states or not 
they are responsible for the screening potential. 

One may suspect that the massive “Schwinger like” states are an artifact of 
the abelianization of the theory in the large Nf limit. To exclude this possibility 
we have found truly non-abelian solutions of the equations of motions. 

The fact that an abelian nature is not necessary for the existense of the massive 
modes tells us that it is plausible that the spectrum is characterized by a smooth 
transition of the screening behaviour from large Nf down to Nf = 1. It may seem 
that there is a contradition between the seminal work of’t Hooft and the various 
types of evidence of the screening behavior of the massless theory. The reason 
for that is that the analysis of ref. [8], which predicts a conhning spectrum, is 
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insensitive to the question of whether the quarks are massless or massive. The way 
to reconsile the two pictures of the massless model is the following. Assume that 
the potential is of the form of eq. (15). In the approach of [8], e^N(. is kept hnite 
in the large limit. This implies that the potential behaves like (1 — e ~ 

-^^\R \[(1 — + o(;^)] for hxed R and large Nc, with k = 2y/2acNc a hnite 

constant. Now it is clear that in the limit of A^c oo the potential looks like a 
linear potential which obviously admits a conhnement behaviour. Thus, the large 
Nc limit prevents one from detecting the truely screening nature of the massless 
system. 

In the present paper we have not discussed the solutions of the equations of 
motion and the corresponding potentials for the case of massive quarks (for large 
Nf it is discussed in [12]). However, using the analogy with the massive Schwinger 
model one can get a general picture of the passage to a conhnig behavior. The 
mass of the massive state of the Schwinger model is shifted once quark mass is 
turned on. But an additional massless state emerges. Exchange of the latter mode 
causes conhnement. Presumably a similar situation occours in the non-abelian 
case. Moreover, in the double bosonization description of the massive model where 
the external sources are also bosonized, the conhning nature of the theory manifest 
itself via the absense of quark soliton solutions. This hts well with the conhning 
spectrum discovered by’t Hooft, the baryonic spectrum analayzed in ref. [7] and 
the results of ref. [2]. One can envisage having a term m?\R\ in the potential, 
where m is the quark mass, in addition to the screening term, with only a screening 
remaining for m = 0. 

Several open questions arise following the results of the present work, for in¬ 
stance the quantization of the non-abelian solutions of the equations of motion, 
deriving similar solutions for the massive model, etc. In the introduction the long 
standing puzzle of conhnement in non-abelian gauge theories in 4D was mentioned 
as one of the motivations for the present work, so we cannot end this discus¬ 
sion without addressing the question of the relevance of our results to the real 
world QCD. The screening of external quarks in the fundamental representation 
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by dynamical fermions in the same representaion both in the massless and massive 
cases seems to £t the pictnre one has abont 4D QCD with fnndamental qnarks. 
The analysis of 2D model with adjoint fermions exhibit a screening behaviour for 
massless quarks and conhnement one for massive case. In 4D one believes that 
adjoint quarks cannot screen external fundamental ones. The mass of the dynam¬ 
ical fermions does not play any role in this issue in 4D. This phenomenon that 
the nature of the mass term in 2D is very different than the one in 4D, was found 
also in other circumstances like the baryonic spectrum. One may speculate that 
the analog of quarks in 4D, whether massive or massless, are massive quarks in 
2D. In [2] it is speculated that the analog of the phase transition that occours by 
turning on mass in the adjoint case is that of breaking SUSY and loosing screening 
in Y = 1 supersymmetric YM in 4D. 

Acknowledgements: One of us (J.S) would like to thank A. Zamolodchikov for 
numerous fruitful discussions. 
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APPENDIX A 

Non-abealian solution of the equation of motion for SU (2) 

Let us parametrize the SU(2) gauge held in the following form A = ff + fr + 
foTQ where /o is real, / is the complex conjugate of / and the matrices r obey 
the algebra [ro,r] = r; [ro,f] = —r; [r, r] = 2tq. In terms of those variable the 
equation of motion eqn ( 8 ) read ( with Jext = 0 ) 

{dd + dc)/ - Afdfo - fodf) = 0 

- - (A.l) 

{dd + ac)fo + 2i{dff - fdf) = 0 

Using the notation / = eqn. (A.l) takes the form 

{dd + dc)g - gd9{fo + d9) = 0 

gdd9 + d9dg + d9dg - {gdfo - fodg) = 0 (A. 2 ) 

{dd + dc)fo + 4:d9g‘^ = 0 

The case where d9 = d9 = 0 is the abelian solution since all the commutator terms 
vanish. In case that d9 = k] d9 = k where k and k are constants it is easy to check 
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that the equation do admit a non-abelian solution as eqn. (12) of the form 


9 = kz + kz] /o 


ar 


kk 


9 


{kk — ac)ac 

4P 


(A.3) 


This result follows also from the ansatz A = e (Aq constant) for which 

eqn. (8) takes the form 


-d9de[To, [ro, Aq]] + idde[To, Aq] + dc^o + d9[[TQ, Aq], Aq] = 0 (A.4) 

inserting the values for the commutators one hnds otcf^+^dOg^ = 0, namely, 89 = k 
where /c is a constant. Using this property leads to a similar relation 89 = k where 
k is also a constant and to the determination of g and /o as given in eqn. (12). 
When an external source of the form Jext = tqJq is added the only modihcation of 
eqn. (A.2) is that the r.h.s of the last eqnation becomes acJo- Similarly, the r.h.s. 
of (A.4) becomes acJext- When Jq is not a constant, bnt a fnnction of f, we need 
to add to 0 a function of 0 , as discussed in section 5. 


APPENDIX B 

Derivation of the Energy-momentnm tensor of QED 2 

In order to nnderstand the energy-momentnm tensor of QCD 2 -, we start with 
the derivation in the abelian theory. The Lagrangian density of massless QED 2 is 
given in bosonized form by 

k^QED = 84>d(t) + \e‘^ + [8(j)A - (B.l) 

where E = Ez^z = 8A — 8A. The corresponding equations of motion are 

28d(t) = — 8F = -^8(f)-, BE = -^d(p (B.2) 

VTT 

We have here the currents j[u{i)] = and i[u{i)] = The definition of 


22 



the energy momentum is 


T - — 

ZZ — 




d{d(j)) 


d{dA) 


= d(j)d(f) + F{dA) H— -=d(i)A — C 

VTT 

= \F'^ + d{FA) 


(B^3) 


In a similar way we have that 


T - — 

■ i - 7 . 7 . - O 


d[FA) 


Tzz = d(j)d(p H— -=d(j)A 

/TV 


Tz 7 = dddd) - 1=8 (f) A 

'TT 


(B.4) 


Note that we have here the “canonical” energy-momentum, not the symmetric one, 
but it generates too the Poincare group. 

In the gauge yl = 0 the above expressisons become 


T - — - F 

■L 77 - rt 1 


1 IT'S. 


Tzz = - d{FA) 


Tzz = d(j)d(j) H— -=d(j)A 

A 

Tzz = d(j)d(j) 


(B.5) 


We now take 2 ; as our “space” variable and 2 ; as “time”. Thus 


P = j dzTzz = j dz{d(j)f = TT J 

P= j dzT^s = / 


(B.6) 


In the non-abelian case, in the gauge A = 0, Tzz has no contribution from 
u~^duA term nor from the term. Thus it is the same as the one derived from 
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the ungauged WZW action, namely 


Nf + Nc''^ ' 

Tzz has no contribution from the WZW part of the action, and the contribution 
from uT^duA is cancelled by the same term in —Thus there is a contribution 
only from term, ie. 

T- — ip2 
-'■22 — 2 "^ 

The corresponding light-cone component P is given by 


P = 



dz 


TT 

Nf + Ne 



For P we need Tzz since P = f dzTzz, but as the difference with the symmetric 
energy-momentum tensor does not contribute to P we also have P = f dzTzz = 
f dz^P^. Note from the explicit expression for the abelian case that indeed Tzz 
differs from by a <9 term that does not contribute to the integral dehning P. 


APPENDIX C 

Minimum of energy for soliton solutions 

The expression of the Hamiltonian for static conhgurations of the (ungauged) 
massive WZW model is 


H 


d‘^zTr[digdig^] + m^Tr[2 — (g + g^) 


(C.l) 


This is also the massive non-linear a model. Consider a diagonal solution ^tq of 
the equation of motion. A general non-diagonal conhguration can be written as 
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g = where B is a U{N) matrix that is time independent. Substitution of g 

into (C.l) one dinds 

H = j dhTr[{digo + [B^diB,go]){digl + [B^diB,gl])+m^Tr[2-{go + gl)] (C.2) 

The difference between the latter Hamiltonian and that which corresponds to ^'o is 

H-Ho = j dhTr[[B^diB,go][B^diB,gl] + 2B^diB{gldigo + godigo])] (C.3) 

For a diagonal solution ^'o fhe last expression reduces to 

H - Ho = j dhTr I [B^diB, go] [B^diB, } > 0 (C.4) 

Thus if g is non-diagonal, it cannot be a classical solution, as after diagonalization 
to go it will have a lower energy. 
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